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In this work, we calculate the matrix elements for the Oz^/3/3-decay of using the deformed 

pn-QRPA method. We adopted the approach introduced by Rodin and Faessler [Phys. Rev. C84, 

014322 (2011)] and by Simkovic et. al. [Phys. Rev. 087,045501(2013)] to restore the isospin 
symmetry by enforcing MJ." = 0. We found that with this restoration, the Fermi matrix elements 
are reduced in the strongly deformed by about 15 to 20%, while the more important Gamow- 

Teller matrix elements remain the same. The results of an enlarged model space are also presented. 

This enlargement increases the total (Fermi plus Gamow-Teller) matrix elements by less than 10%. 

PACS numbers: 23.40.Hc,23.40.Bw,27.70.+q 


The neutrinoless double beta decay (hereafter QvjSfi- 
decay), if it exists, is one of the rarest processes in our 
universe. Its discovery could lead to the dawn of a era 
for new physics beyond the Standard Model. Various ex¬ 
periments search or are proposed to look for this exotic 
process. Currently, about a dozen isotopes have been 
confirmed to decay by the two-neutrino double beta de¬ 
cay (2z//3/3) mode (for a review see i)- They are also 
good candidates for Oi//3/3-decay. Of these isotopes ^^°Nd 
has a relative small 2vfif3 half-life as shown in Table I of 
ref. [I|. This decay has also a large Q value (Q = 3.367 
MeV, see table I of ref. i) . Recent calculations in ref. 
3 show, that it also has the largest phase space for the 
light Majorana neutrino mechanism - the most proba¬ 
ble mechanism for this decay. Thus this decay system is 
a promising candidate with possibly the shortest 0^/3/3- 
decay half-life. To confirm this assumption, one needs 
to investigate further this process. Compared to other 
/I/3-active isotopes, ^^°Nd is supposed to be strongly de¬ 
formed. This poses difficulties for exact shell model cal¬ 
culations. So for this nucleus one needs to resort to other 
methods, e. g.: the Projected HFB [^, the Interacting 
Boson Model (IBM) Q and Energy Density Functional 
(EDF) (non-relativistic ref. @ and relativistic ref. 0)- 
These methods calculate the ground states of the initial 
and final nuclei of the double beta decay, and then use 
closure for the intermediate states with an averaged en¬ 
ergy denominator to obtain the nuclear matrix element. 
In this way one needs not to calculate the wave functions 
of the intermediate nucleus ( ^®°Pm). Another category 
of methods used extensively is the Quasi-particle Ran¬ 
dom Phase Approximation (QRPA), which calculates ex¬ 
plicitly the intermediate states. According to the mean 
fields and interactions, different versions of QRPA are 
in use. Our group uses QRPA with realistic forces [l^ . 
There exist also QRPA calculations with the Skyrme in¬ 
teraction [ 1 ^ and its variations m- 


QRPA calculations. With this new formalism, isospin 
conservation is obtained by setting the value of Mp'' to 
zero. The results in ref. 0 show, that for spherical 
nuclei this new approach reduces the Qv(3j3 Fermi ma¬ 
trix elements by about 30^40%, while the more impor¬ 
tant Gamow-Teller contributions are unchanged. In this 
work we adopt this new formalism for the ,5/3-decay of 
strongly deformed ^®°Nd. This is done by separating the 
renormalization of the particle-particle residual proton- 
neutron interaction into the T=I {Opp^) and T=0 
parts. The new treatment differs from our work in [I2j| . 
where Qp^^ = Jp^^ and the gpp's are fitted to the exper¬ 
imental two neutrino double beta decay matrix elements 
M'qp. The previous approach [T^ yields a relative large 
value of the Fermi part Mjf', which should disappear 
according to isospin conservation. In the present work 
we restore (at least partially) the isospin symmetry for 
the first time in the QRPA approach for deformed nu¬ 
clei and calculate the Oiz/3,5-decay matrix element for the 
deformed ^®°Nd nucleus. 

At first we give a brief review of our method. The 
QRPA states are defined as: 

\K\ m) = J2 - YpZanap) |0) (1) 

pn 

(Bp Byj = A"), '^p'^n — 


Where the a’s are the quasi-particle creation and anni¬ 
hilation operators. is the angular momentum projec¬ 
tion to the symmetry axis of the axially symmetric de¬ 
formed nucleus. |0) is the QRPA vacuum. But in actual 
calculations it is always simplified to the BCS vacuum. 
X’s and Y’s are derived from the solutions of the QRPA 
equations in the deformed system as in ref. 13|: 


( A{K) B{K) 
{-B{K) -A(K) 



( 2 ) 


Recently, a new formalism ii with restored isospin 
symmetry has been developed by the Tubingen group for 


Now due to partial isospin restoration the detailed ex¬ 
pressions of A and B for the realistic G-matrix are 
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slightly different from those in 
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Here the particle-hole interactions Vpnp'n' are the same as 
in [T^ . The particle-particle interactions are now isospin 
dependent. They are now different for the T=0 and T=1 
parts. These two parts are expressed by expansion into 
the G-matrix elements in the spherical basis: 

yT=o{i) Y - 

pnp'n' / ^ ^ PVp''^Vn p'Pp/n'p^/ 

J VpVriVp'Vn' 

X G'^^°^i{'qpgnVp’gn'J) (4) 

The decomposition coefficients F are defined in [l^ . 

QRPA yields the matrix elements of the 2vf3P- and the 
Oi^/3/3-decay as: 


= E {pnK^Wn'K^)raj,m. 

mi ,m f 

X {0'^\clcn\K'^rnf){K^mf\K'^rni) 

{K^m,\clcn\0t) (5) 

For the 2j//3/3-decay, the hrst term in the above formula 
can be divided in two parts: 


{pnK^\Ol^T\pnK^)n.,,m. 

{pnK^\Of\p'n'K^),^^,m. 


(p|r“''cr|n) • (p'lr+crln') 
Eg.s. + {Enii + 

(p|r+|n){p'|T+|n') 

Eg.s. + {Errii + Emf)/^ 

( 6 ) 


The detailed expressions of the single particle transition 
matrix elements (plr ^gln) in the deformed system can 
be found in [l^ . Il3j |. Similar expressions can be ob¬ 
tained for the F matrices. The denominators are now 
slightly different from [l^. Here Eg ^, = (2M(^^°Pm) — 
M(i50Nd)-M(i50Sm))/2, and Fl^’s are given in ref. jl3 
as Em = ujm — Wist- The w’s are eigenvalues of QRPA 
and uJist are the lowest eigenvalues of QRPA for different 
AT’^’s. For QRPA, the transition matrix elements between 
the ground states and intermediate states are expressed 
as: 


(A'’",to|cJ,c„|0+) = {UpVnXp^"^ - VpUnY^'^) (7) 

m’s and v's are the occupation amplitudes of the BCS 
equations. X’s and Y’s are solutions of the QRPA equa¬ 
tions. The overlap of the initial and final intermediate 
states in eq. ([S]) is more complicated and is given in refs. 


TABLE I: A brief summary of the parameters used in the 
current calculations which are changed compared to ref. [T^ . 
Other parameters are not changed. The last line displays 
the old parameters used in ref. [T^], where both gpp values 
for T = 0 and T = 1 are the same and thus isospin is not 
restored. 



A = 4 - 6 

QA = 1.27 QA = 1.0 

N = 0-7 

QA = 1.27 PA = 1.0 

Qpp 

Qvv^ 

1.08 1.03 

1.34 

0.77 0.74 

1.06 

ppp[12] 

1.05 1.00 



[3 EE [ 3 . The overlap between the two BCS vacua 
f{BCS\BCS)i suppresses the matrix elements strongly, 
if the deformations of the initial and final nuclei are differ¬ 
ent lig. W e discussed this effect on the matrix element 
in ref. |l7j| . In this work we take the above BCS overlap 
equal to 0.52 as given in Table I of ref. El. 

The 0i^/3/3-decay operators are two body operators with 
integrations over the loop momentum q, see [iE.EE|. The 
total Ovfjp matrix element can be expressed as: 

It contains a Fermi (F) and a Gamow Teller (GT) part 
and depends on the vector gv and the axial vector gA 
coupling constants. M^'^’s contain summations over dif¬ 
ferent transition operators EE- In this work we neglect 
the tensor contributions. Th^ are usually small com¬ 
pared to the other two parts Q. 

We now give briefly the parameters of the model used 
in this work. The single particle energies and wave func¬ 
tions are obtained by solving the Schrodinger equation 
with a Coulomb corrected Woods-Saxon potential. The 
same deformation parameters (32 , for ^®°Nd (32 = 0.240 
and ^^°Sm (32 = 0.153, as in ref. El are used. The pair¬ 
ing strengths gpaiVs are obtained by fitting the experi¬ 
mental gaps. The strength of the particle-hole residual 
interaction gph is the same as in our earlier description 
of the Ov(3(3 decay in ^®°Nd. The gpp now changes in the 
T=1 channel see table [H They are obtained by fit¬ 

ting the experimental 2^/3/3-decay matrix elements. We 
found a similar behavior as in ref. Q : The GT part Mq^ 
depends only on gpj("° and not on While the Mp’' 

depends strongly on g^^^ but not on gp(p^. Thus gpp^ 
and Ppji^ can be determined separately by reproducing 
the experimental GT matrix elements, which for these 
nuclei are = 0.07 [l| and by enforcing Mp'^ = 0 as 
required by isospin symmetry. (The reason that we call 
this ’’partial restoration of isospin symmetry” is, that in 
exact shell model calculations, where the isospin symme¬ 
try is conserved, all the successive single particle transi¬ 
tions from initial to intermediate and then to final states 
are exactly zero. For QRPA only the overall sum of these 
transitions disappears.) Ref.jE shows, that g^""^ should 
be approximately the same as the pairing strength dpp 
and dnn in the T = 1 channel. We checked this in the 
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present calculations in the large model space. We find 
for ^®°Nd, dpp = 0.94 and d„„ = 1.03, and for 
dpp = 0.95 and dnn = 1.04^ This gives an average pairing 
strength of d = 0.98 and d = 1.00 respectively, about 5% 
smaller than gp^^. These results agree with ref.@ and 
imply that the new parameter introduced is consistent to 
the pairing strength in the T=1 channel. 

In the present calculations, we fit two sets of val¬ 
ues, one for the bare = 1-27 and another for the 
quenched gA = 0-75gAo values obtained from experiment 
(l^ . Using calculated 1^ excitation energies in the en¬ 
ergy denominator, barely changes the matrix elements 
compared to ref. [T^ . The newly fitted gp^^ values only 
differ by 1 to 2 % while M^j. is basically not changed. 
As a result of improved computational facilities we can 
now use a much larger model space of up to eight major 
shells N = 0 — 7 compared with a model space = 4 — 6 
used in ref. [H- The values of gpp in this larger model 
space are smaller. This implies that in a smaller model 
space, the interaction in the particle-particle channel is 
enhanced to account for the smaller model space. The 
enlargement of the model space changes the results of 
0z//I/3-decay as we shall show below. 

We illustrate the results of single intermediate K'^ con¬ 
tributions in FiglT] The blue bars are the results in the 
small basis N = 4 to 6 without isospin symmetry restora¬ 
tion and with the conventional choice of 

^PP ^ ^PP^ 

The symmetry restored results in the small basis are dis¬ 
played by the red bars. Here as for all results calculated 
in this paper the CD-Bonn nucleon-nucleon Brueckner 
G-matrix elements and the Brueckner short range corre¬ 
lations (src) of ref. are used. 

The effect of the isospin restoration leads to changes 
of Mp'' for each K^. The largest change is obtained for 
= 0+, which corresponds to J'^ = O’*' in the spher¬ 
ical limit in ref. For K'^ = 0+ the Fermi matrix 

element is reduced to about 1/3. Changes for other 
AT’^’s are not that significant. Compared with the con¬ 
ventional parametrization, the new formalism increases 
by 0.35 (more than 30%, see table. HI. But this large 
increase does barely change the values of Mq'p, since its 
main contribution is due to T = 0 nucleon pairs. Only 
Mp' for = O'*' is sensitive to gp^^, because the main 
contribution originates from the interaction in T = 1 nu¬ 
cleon pairs. 

In fig. [TJ we show also how the enlargement of the 
model space affects the final results. Amazingly pre¬ 
vious truncations of the model space, though numeri¬ 
cally insufficient due to its small size, produce however 
similar results as the ones obtained from a large model 
space. The main increase of M^p is due to the two states 
— O'*' and = I’*' by the larger model space. For 
other iF^’s the large model space could either increase 
or reduce slightly the matrix elements depending on the 
detailed transitions. This doesn’t mean that the contri¬ 
butions from transitions outside of the truncated model 
space N = 4 to 6 are not important, since we have quite 
different renormalization strength parameters 5 pp’s for 


TABLE II: A summary of the results calculated with the 
Bonn CD potential with different model spaces, with or with¬ 
out Brueckner short range correlations (scr and no scr) of 
the Bonn CD nucleon-nucleon interaction, two axial vector 
coupling constants qa = 1.27 and qa ~ 0.95 and with and 
without partial isospin restoration (ISR, no ISR). The matrix 
element M'^’' is defined in eq. (IS)). 



gA = 1.27 

Mg'/ M'°" 

qa = 0.95 

Mg'^ Mg/ M'°'' 

N=4-6 ISR, no src 
no ISR, no src 
ISR, src 
no ISR, src 

-1.308 2.081 2.891 
-1.565 2.091 3.061 
-1.367 2.214 3.062 
-1.628 2.224 3.233 

-1.306 2.371 2.143 
-1.614 2.381 2.340 
-1.365 2.508 2.257 
-1.679 2.518 2.457 

N=0-7 ISR, no src 
ISR, src 

-1.390 2.309 3.171 
-1.454 2.466 3.367 

-1.369 2.629 2.328 
-1.433 2.790 2.458 


both T=0 and T=1 parts for the truncated and large 
model space (see table |I|. For the 2 j^/ 3/3-decay larger val¬ 
ues for gpp^s are required for the small model space N = 4 
to 6, because some of the correlations are missing in the 
truncated model space and we need to compensate these 
correlations in the QRPA calculations by increasing the 
interaction strength. The situation for the Oz^/3/3-decay 
is similar to the 2i/ decay. The larger gpp^s in the small 
space mimic the behavior in a much larger model space 
and produce values close to the values of the matrix el¬ 
ements in the large space. So the larger force strength 
gpp compensates for the smaller model space and finally 
one obtains very similar transition matrix elements for 
the Oi//3/3-decay (see table H] and ref. [l^). In general 
(See fig. I) the larger model space enhances the tran¬ 
sition matrix elements calculated with the short range 
Brueckner correlations using the CD-Bonn force (src, see 
also table m especially for GT matrix elements. For the 
Fermi (F) part the situation is different. This indicates a 
smaller sensitivity of Fermi matrix elements on the size 
of the model space. 

In table El we have summarized all the results with 
different model spaces and compared the new formalism 
with Isospin Restoration (ISR) with results without ISR 
in the small model space N = 4 — 6. The isospin restora¬ 
tion (ISR) reduces the Fermi marix elements Mp' by 
about 15 to 20%. The main reduction originates from 
= O'*". M/jy is almost not changed (see table El box 
N = 4-6). The overall matrix element M'^'^ is reduced 
by less than 10% including the ISR. For the dependence 
on the axial charge gA the Gamow-Teller matrix element 
M^P increases for a decreasing gA- The result for the 
total matrix element going from gA = 1.27 to 0.95 is a 
reduction by about 30%. The effect of short range cor¬ 
relations (src) is the same as in ref. [l^. The Brueckner 
short range correlation src for the Bonn CD force slightly 
enhance in the small and the large model space the ma¬ 
trix elements by about 5%. The CD-Bonn interaction is 
known to yield small short range correlation (src). Thus 
one expects larger changes for other short range correla¬ 
tions. 







4 



FIG. 1: (Color online) The decomposition of Qvf5j3 matrix elements for different values of . Here, ”Msp” is an abbreviation 
for "model space”, where "Msp I” refers to the N=4-6 and "Msp 11” to the N=0-7 space. ”ISR” abbreviates "Isospin Symmetry 
Restoration”. 


TABLE III: Results for the total matrix element M'°'^ (defined in eq. (dl). qa = 1.27 and the self-consistent Brueckner CD- 
Bonn short range correlations 0 are used. One should be aware different methods use different conventions such as nuclear 
radii and short range correlations, etc, which affect the final results. Here ’’Non-closure” means, that the intermediate states are 
calculated explicitly, and ’’Closure” means, one calculates the transitions from ground states to ground states without taking 
into account explicitly the intermediate states. 


Non-closure 

previous [12] this work QRPA-SKflT] LP-QRPA[lll 

Closure 

PHFB[^ IBM-2 NREDF[6] REDF[7] 

3.34 3.37 2.71 3.60 

3.24 2.67 1.71 5.60 


The results for the absolute values of the Fermi and 
the Gamow-Teller (GT) matrix elements are enhanced 
by about 5 and 11%, respectively, enlarging the model 
space for both axial charges considered. This enhance¬ 
ment, as we have seen from fig. [TJ stems mostly from 
the K'" = O’*' and K'" = I’*' (J'"’ = I’*') states for the 
GT part. This is due to the sensitivity on for these 
K'^’s. A detailed analysis of these sensitivity will be pre¬ 
sented in a future work. In refs. one argues, that 
in a spherical system this enhancement will be partially 
compensated by a decrease of the Tensor part, which is 
not included in the current calculation. The spherical 
calculations in ref. [2l| suggests, that the enhancement 
due to the enlarged model space could be smaller than 
shown here. Further investigations are needed for such a 
conclusion in a deformed system. 

In table uni we compare recent results of 0j^,5/3-decay 


for from different methods. Gompared with our 

previous results, we now have partially restored isospin 
symmetry, which reduces Mp^. Due to the new energy 
denominator and the enlarged model space, the values 
are slightly increased compared to QRPA calculations 
based on Skyrme forces [0], which do not include short 
range correlations. Our results are smaller than that of 
ref. [ll| , which uses pp and nn QRPA. Except for QRPA 
most methods use the closure approximation, where one 
does not need to determine the intermediate states, but 
treat it by closure. It has been suggested by that 
this approximation changes the results by at most 10% 
for shell model calculations. Some further investigations 
are needed to see, if this holds for all methods. In table 
EH we have listed results from closure methods. They de¬ 
viate strongly from each other. The PHFB gives a final 
matrix element very close to ours. IBM-2 gives a much 
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smaller value close to Skyrme QRPA. The two results us¬ 
ing the enerey density functional (EDF) give the largest 
(5.6, ref. [3]) and the smallest values ( 1.71, ref. p). 
They are different by a factor 3. These authors need to 
investigate further the strange large difference between 
the two ’’density functional” results. 


In conclusion, currently we still have variations of the 
0://3/3-decay matrix elements for different methods by a 
factor of three or more, which can produce an uncertainty 
of the order of one magnitude for the half-lives. Further 
comparisons among methods should be studied to find 
the reason for such large deviations. 
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